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Abstract 

Dynamic properties of a one-dimensional probabilistic cellular automaton are studied 
by monte-carlo simulation near a critical point which marks a second-order phase tran- 
sition from a active state to a effectively unique absorbing state. Values obtained for 
the dynamic critical exponents indicate that the transition belongs to the universality 
class of directed percolation. Finally the model is compared with a previously studied 
one to show that a difference in the nature of the absorbing states places them in 
different universality classes. 
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1 Introduction 



Discrete models of nonequilibrium stochastic processes form a class of interacting particle 
systems [|l|. Of the models studied with short range and translationally invariant interactions 
in space and time, the ones exhibiting a continuous phase transition from an active steady 
state to an absorbing state fall into two universality classes : 

(1) the class of directed percolation (DP), 

(2) the class of parity conservation (PC). 

Models with a unique absorbing state have been conjectured to belong to the DP class p, 
^. This is yet the larger of the two classes and includes, for example, lattice models of di- 
rected percolation in d + 1 dimensions 0, the contact process for an epidemic 0, Schlogl's 
first and second models of autocatalytic reactions 0], the Domany-Kinzel automaton ||^, |^, a 
lattice version of reggeon field theory [TO] and branching annihilating random walks with an 
odd number of offsprings ||ri|. The order parameter is usually the density of 'particles' (oc- 
cupied sites, kinks). The models mentioned above are one-component systems and therefore 
the order parameter is scalar, a requirement of the DP-conjecture P, §]. The DP-conjecture 
was generalised to include multicomponent systems such as the ZGB model of heterogeneous 
catalysis [jl2|, |13|, and interacting dimer-trimer models The class of DP was shown 
to further include systems with an infinite number of absorbing states these absorbing 
states are frozen configurations characterised by unique statistical properties, a lack of long 
range correlations and a general lack of symmetry among them. 

On the other hand, models with degenerate, mutually symmetric absorbing states are 
believed to belong to the parity-conserving class where the number of 'particles' are conserved 
modulo 2. Prominent examples are the probabilistic cellular automaton models A and B 
of [T^, |TB|, an interacting monomer- dimer model and branching annihilating random 



walks (BAW) with an even number of offsprings |2y] . According to one point of view ||2|, |T 
the mechanism that puts these models in a class different from that of DP is the conservation 
of particle number modulo 2. This point was proved for BAW with an even number of 
offsprings: introduction of spontaneous annihilation of particles in the model destroyed the 
conservation of their number modulo 2 and the critical behaviour of this modified model 



was observed to be in the class of DP ||2l[ . According to a second point of view the critical 
behaviour of the models in the PC class is due to the symmetry among its absorbing states 



22, 23|. To prove the point the interacting monomer-dimer model of [19|, in the presence of 
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a weak parity-conserving field that destroyed the symmetry among the absorbing states, was 
shown to exhibit critical behaviour in the DP class . This view was further emphasized 
in the generalised versions of the Domany-Kinzel automaton and the contact process [p^] . 
In these generalised models there was no explicit parity-conservation law and with two 
symmetric absorbing states the critical behaviour was in the PC class. In the presence of 
a symmetry breaking field the critical behaviour of these models changed to the DP class. 



However in the case of the probabilistic cellular automata of |Tj, [T8[ both the conservation 
of particle number modulo 2 and mutual symmetry among the absorbing states are present 
and therefore it it is not apparent which of these two features is responsible for the models 
to be in the PC class. 

In this paper I shall study, using time-dependent monte-carlo simulations, the dynamic 
critical properties of a one-dimensional probabilistic cellular automaton which has three 
absorbing states and exhibits a phase transition from a active state to one of them only. The 
simulations provide values for the critical point (more accurate than previous estimates |^5|, 
p6| ) and the dynamic critical exponents that decide the universality class to which the phase 
transition belongs. 



2 The model 

The probabilistic cellular automaton studied here is 'elementary' (in the sense of Wol- 
fram [|^) with two states per site and translationally invariant nearest neighbour inter- 
actions. The probabilistic behaviour enters the model through two mutually symmetric 
components of the evolution rule (like specific noise added to a otherwise deterministic sys- 
tem), the rest of the rule components being deterministic in nature. 

Formally the model is defined as a line of sites with a binary variable Xi G {0, 1} 
assigned to each site i. A site is said to be occupied if = 1 and unoccupied otherwise. 
Starting from a given configuration {a;^'''*} the system evolves by parallel update of the 
variable Xi at all lattice sites following a local rule of evolution. With nearest neighbour 
interactions the evolution rule is defined by a set of eight components [xfli, xf\ xf^i] t— >• 
xf^^'' corresponding to 2^ distinct three-site neighbourhoods : 
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1 with probability p 
with probabihty 1 — p 



The evolution rule (|l|) thus follows, according to Wolfram's nomenclature scheme 'ele- 
mentary' rule 4 with probability 1 — p and 'elementary' rule 22 with probability p. 

Clearly, the components of the evolution rule are of two kinds: (1) active components, 
where the central site changes its value, and (2) passive components, where the value of the 
central site remains unchanged. The dynamic evolution of the system is due to the active 
components. This involves two opposing processes : 

(a) Annihilation of adjacent occupied sites (multi-'particle' annihilation) due to the rule 
components 111h-^0, llOi-^O, Ollt-^O prevents the survival of occupied pairs of nearest 
neighbours; this is a deterministic process. 

(6) Creation of an occupied site (100 i-^ 1, 001 h-^ 1) requires an unoccupied site to have 
exactly one occupied neighbour; this process of creation occurs only with a probability p. 
The multi-'particle' annihilation and creation of 'particles' can, in effect, give rise to a 
diffusion process — if a 'particle' at site i gets annihilated after creating another 'particle' at 
a neighbouring site i + r, it has effectively taken a step from i to i + r. 

The passive components determine the absorbing states of the model. Since there is 
no spontaneous creation of occupied sites (000 ^— 0) the 'vacuum' (all sites unoccupied) is 
always an absorbing state : 



Absorbing state I : = for all i. (2) 

Again, it is evident from the rule components 010 ^ 1 and 101 i— > that an occupied site 
with unoccupied neighbours remains occupied and vice versa. These features of the evolution 
rule lead to two mutually symmetric absorbing states : 

Absorbing State II : Xi 
Absorbing State III : Xi 



for i = even, , . 

1 for i = odd. ^ ' 

1 ioY i = even, , . 

for i = odd. 



The main point of interest in this paper is a phase transition exhibited by the model [P5|]. 
For p less than a critical value pc there exists three distinct steady states given by the three 
absorbing states of the model. For all initial states but two, it has been observed in computer 



simulations that the only steady state is the 'vacuum' (absorbing state I). The two cases 
of exception occur when the initial state is either absorbing state II or absorbing state III, 
which must also be respectively the steady states of the system. Evolving from any other 
initial state the mutually symmetric absorbing states II and III are never reached. This 
can be understood by the fact that these states are fixed points that repel — a damage 
introduced in these two states by flipping only a single bit spreads through the entire lattice 
and eventually the 'vacuum' is reached. Above pc there is another steady state called the 
'active' state with a constant non-zero density of occupied sites. The active state is stable 
only on an infinite lattice; for finite lattices it occurs as a metastable state that will decay 
to the 'vacuum' if allowed to evole for sufficiently long time. On an infinite lattice, in the 
supercritical region [p > pc), all possible initial states other than the three absorbing states 
evolve into the 'active' steady state. The density of occupied sites poo in the steady state 
acts as the order parameter for this phase transition — in the supercritical region poo goes 
to zero continuously as p approaches pc pSf : 



Poo^{p-Pcf, P^Pc + , (5) 

where (3 is the critical exponent for the order parameter. For random initial states the model 
therefore makes a continuous (second order) phase transition, at p = pc, from a active state 
with Poo > to a effectively unique absorbing state, the 'vacuum'. Because of the conjecture 
of P, this phase transition is expected to belong to the universality class of DP. 

3 Dynamic properties at the critical point 

The phase transition is characterised here by critical exponents describing the dynamic 
properties of the model at the point of transition (the critical point pc). To do so dynamic 
properties of the model are studied by monte-carlo simulations on a computer, only close to 
Pc- While the study of steady state properties require simulations starting from disordered 
states, the dynamic properties are studied using initial states with a single occupied site. 
Following the evolution rule (^ an initial occupied site grows into a cluster; the position of 
this initial occupied site is called the 'origin' of the cluster. For each value of p, 10^ clusters 
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were simulated. Each cluster was allowed to evolve for 5000 time steps, unless it had died 
out earlier. Typical examples of evolution near the critical point are shown in figure 1. 

The quantities measured are : (1) the survival probability P{t), which gives the chance 
that there is at least one occupied site after t time steps, (2) the average number of occupied 
sites N(t) after t time steps, and (3) the mean square radius -R^(t) of the cluster (or, the 
mean square displacement from the 'origin' of the cluster) after t time steps. At the critical 
point p = pc, these quantities are expected to follow power-type scaling laws in the long-time 
limit (t ^ oo) : 



Pit) ~ r', 

N{t) ~ t\ (6) 

where 6, rj and z are dynamic critical exponents. In the case of N{t) the average is taken over 
all clusters including those which have died out, while -R^(t) is averaged over the occupied 
sites in the surviving clusters only. 

Results for the three quantities P(t), N(t) and R'^it), obtained from computer simulations 
of the model (|lj), are shown in figure 2. On log-log plot curves in the subcritical region bend 
downward while those in the supercritical region bend upward; at the critical point the curves 
are expected, if the scaling laws are true, to be straight lines as t ^ oo. It is obvious from 
figure 2 that 0.75 < Pc < 0.753. More precise estimates for pc and the critical exponents are 
obtained by the method of effective exponents due to Grassberger Effective exponents 



6t, rjt and are defined as the local slopes of the curves shown in figure 2 : 

_ ^ A[logP(t)] ^ A[logiV(t)] ^ A[\ogR\t)] 

* A[logt] ' A[logt] ' * A[logt] ' 

which are measured by using the formulae 



_ ^ log[P(t)/P(t/m)] ^ 
logm ' 

and similar expressions for rjt and Zt. Results for m = 10 are shown in figure 3. Like the 
curves in figure 2, curves for p < p^ bend downward while those for p > p^ bend upward. 
For finite times the dynamic quantities do not have a pure power-law scaling form like 
in general there exists correction-to-scaling of the type : 
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P(t)~r^(l + ^ + ^ + ---), (9) 

and similar expressions for N{t) and -R^(t) with correction-to-scaling exponents 6', rj' and z' 
respectively. Consequently the behaviour of the effective exponents defined in (|^) are given 
by0: 

a S'b , ^ 

St = S+- + ^ + ---, (10) 

and similar expressions for rjt and Zt- The critical exponents of the model appear as the 
asymptotic values of the corresponding effective exponents : 



6 = hm^^oo St , 11 = lim^^oo , z = limj^oo Zt . (11) 

Therefore in a plot of an effective exponent (as the ordinate) versus 1/t (as the abscissa) the 
corresponding critical exponent is obtained as the intercept of the curve for p = pc on the 
ordinate axis. Using this method the following estimates for the critical characteristics of 
the model were obtained from computer simulation data : 

= 0.7515 ±0.0005 (12) 

and 



5 = 0.16 ±0.01, 

T] = 0.32 ±0.02, (13) 
z = 1.27 ±0.01. 



The value of the critical point agrees closely with previous estimates |]25|, and improves 
upon them. The values of the critical exponents are also found to satisfy the scaling rela- 
tion [Tg] : 



dz = 2r] + 4:6, (14) 

where d is the number of spatial dimensions of the system (here d = 1). 

Evolving from disordered initial states, the density of occupied sites pt was observed to 



follow the same dynamic scaling as the survival probability. In computer simulations of [£5 
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the exponent characterising the critically slow relaxation pt ~ t"" aX p = pc was found to be 
a 0.16 ^ 5. 

The values of all the three dynamic critical exponents agree, within the limits of error, 
with the corresponding values for DP in 1 + 1 dimensions [l^, In a previous work p6| 



the critical exponents v± and z/y for the correlation length and correlation time respectively 
at this particular phase transition, obtained by finite-size scaling methods, were also found to 
belong to the DP class. However, the order parameter exponent (3 was observed to disagree 
with the DP value ||2^; this was an error arising out of finite-size effects and fluctuations 
due to the small sample size used for averaging. Since 5 and v\\ are already in the class of 



DP, the exponent /3, by virtue of the relation 13 = d^w ||T0[, must also agree with the directed 
percolation value for 1 + 1 dimensions. 

The last result is concerned with the fractal dimension of the clusters in the single space 
dimension at p = Pc- The average number of occupied sites per surviving cluster is given by 
Ns{t) = N{t) / P{t). The fractal dimension dp of the clusters at fixed time is defined by : 

Ns^R'^^. (15) 

Following the definition of the dynamic exponents the fractal dimension is expected to 
satisfy the relation : 

dFZ = 2{ri + S). (16) 
Figure 4 shows a log-log plot of Ng versus R at p = pc- The slope of the curve is given by : 

rfi. = 0.74 ±0.02, (17) 
which satisfies relation (|T6|) within the limits of error. 

4 Discussion 

In this concluding section I shall compare the model defined by ([l|) with another that has 
the same set of absorbing states and that belongs to a different university class. 

The probabilistic cellular automaton (|I]) studied here was found to have three absorbing 
states given by (^, (^) and (^. For p < Pc the 'vacuum' (absorbing state I) is the only 
attractor of the model while the other two absorbing states (II and III) are never reached 
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from disordered initial states. Consequently the 'vacuum' appears, in effect, to be the unique 
absorbing state in the subcritical region. The dynamic critical exponents characterising the 
phase transition in this model indicate that the transition belongs to the class of DP, in 
agreement with the conjecture of |0, 

On the other hand, the phase transition occuring in model A of [|T^ belongs to the PC 



class. This model is yet another one-dimensional elementary probabilistic cellular automaton 
defined by the evolution rule : 



t: 111 110 Oil 101 010 100 001 000 , , 

(18) 



t + 1 : 1110 

^ ^ ^ 

with probability p 

1 with probability 1 — p 



It is remarkable that the three absorbing states of this model are exactly the same as those of 
model (|l|). However, contrary to their nature in model (P, the mutually symmetric absorbing 
states II and III occur as attractors of this model for p < Pc while aborbing state I (the 
'vacuum') is never reached from disordered initial states. It appears that the contrast in the 
nature of the absorbing states between the two models places them in different universality 
classes. In that case the non-DP behaviour of model (|I8]) must be due to the degeneracy in 



the absorbing state in the subcritical region, thus supporting the view of [22, 23 
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Figure Captions 



Figure 1. Typical examples of evolution from a single occupied site following rule : (a) a 
case in the subcritical region {p = 0.74) and (b) a case in the supercritical region {p = 0.76). 
The first 200 time-steps if the evolution process are shown. 

Figure 2. Results for three dynamic properties from monte-carlo simulations of the model 
: (a) survival probabihty, (b) average number of occupied sites and (c) mean square radius 
of the evolving cluster. Each of the three panels contains five curves that correspond to p = 
0.753 (top), 0.752, 0.7515, 0.751, and 0.75 (bottom) respectively. 

Figure 3. The effective exponents measured as the local slopes of the curves shown in 
figure 2. For large t, data have been averaged over many time-steps in order to suppress 
fluctuations. 

Figure 4. The log-log plot of the average number of occupied sites as a function of the 
root-mean-square radius of the clusters at p = pc- The slope of the curve gives the fractal 
dimension of the clusters at fixed time. The straight line drawn below the curve is the graph 
of Ns = const. 
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Figure 1(a) 




Figure 1(b) 
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Figure 2(b) 




Figure 2(c) 
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Figure 3(b) 
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